Superconductivity in the SU(N) Anderson Lattice at f/ = oo 
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We present a mean-field study of superconductivity in a generalized A''-channel cubic Anderson 
lattice at [/ = oo taking into account the effect of a nearest-neighbor attraction J. The condition U — 
oo is implemented within the slave-boson formalism considering the slave bosons to be condensed. 
We consider the /-level occupancy ranging from the mixed valence regime to the Kondo limit and 
study the dependence of the critical temperature on the various model parameters for each of three 
possible Cooper pairing symmetries (extended s, d-wave and p-wave pairing) and find interesting 
crossovers. It is found that the d— and p— wave order parameters have, in general, very similar 
critical temperatures. The extended s-wave pairing seems to be relatively more stable for electronic 
densities per channel close to one and for large values of the superconducting interaction J. 
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I. INTRODUCTION 

The superconducting behavior of heavy-fermion ma- 
terials has attracted rriuch attention due to its non- 
conventional properties .El Despite the large amount of 
work trying to understand heavy-fermion superconduc- 
tivity, the normal state properties, the symmetry of the 
order parameter, the origin of superconductivity and the 
interplay between superconductivity and magnetism are 
still interesting and open questions. 

Some of these materials, such as UAgCu^, UCut, 
U2ZniY, order antiferromagnetically at low temperatures 
while others (such as UBeis, CeCu2Si2, UPt^) order in a 
superconducting state and others show no ordering (such 
as CeAls, UAuPt4, CeCug, UAl2)a There are materials 
which order both antiferromagnetically and become su- 
perconducting as the temperature drops {e.g. URu2Si2, 
C/q. 97^/10.03-6613) and it has recently been found that 
UPd2Al3 shows coexistence of superconductivity and 
local moment antiferromagnetism.H All these materials 
have very large specific heat coefficients 7, indicating 
very large effective masses, hence the designation heavy 
fermions. 

The superconducting properties of a system depend 
on the type of ground state that the system exhibits in 
the normal phase. The large specific heat 7— coefficient 
can have|-.two very different origins: a Kondo-impurity 
behavior Jj in which case 7 behaves as the inverse of the 
Kondo temperature Tr-; or a Kondo-lattice behavior, in 
which case 7 is controlled by a large density of states at 



the Fermi energy. The large density of states arises from 
a hybridization mechanism between the conduction band 
and localized electronic states (/—states, say).Q 

Even though the large effective masses indicate strong 
correlations between the electrons, these behave in many 
cases as essentially "free" , with rcnormalized parameters, 
as explained by the Fermi liquid theory. However, there 
has recently been growing evidence that other material^ 
have properties that do not fit the Fjsrmi liquid picture.El 
The reason could be either disorder J3 vicinity to a quan- 
tum phase transition,l3 or unusual impurity-like behavior 
such as the one described, by generalized models, as the 
n-channel Kondo model.El The n-channel Kondo lattice 
shows interesting behavior and it has been shown to be 
an incoherent metal at low temperatures with a resid- 
ual entropy that is usually lifted via ordering at very low 
temperatures 

A consistent description of the overall properties of the 
heavy-fermion behavior has been achieved assuming that 
a generalization of tiijG-|impurity Anderson model to the 
lattice case is valid£l't2l In the Anderson lattice the en- 
ergy of a single electron in an /—orbital (e. g. 4/^) is 
eo, and the energy of two electrons in the same /—orbital 
(4/-^) is 2eo + U, where U is the on-site Coulomb repul- 
sion. The energy of the 4/^ state is much larger than 
the energy of the if^ state. Moreover, these systems are 



often characterized by 
the spin-orbit coupling 



e angular momentum, due to 
In general, both the large val- 
ues of U and the large total angular momentum must be 
included in any model used to describe the properties of 
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heavy-fermion materials. 

The SU(N) Anderson lattice Hamiltonian is believed 
to give a good description of the normal state of Kondo- 
lattice systems.cl The limit [/ = oo is considered in many 
calculations since the experimental U values are large. 
The Anderson lattice model predicts Fermi-liquid like 
properties in the normal non-magnetic state. The theo- 
retical results give a good description of many materials 
and explain the main features at low temperatures such 
as universality, large effective masses, the Kondo reso- 
nance at the Fermi level. At the single-impurity level the 
picture is clear. In the Kondo limit the /-level has an 
occupation close to one leading to a localized spin that is 
shielded by a conduction electron spin cloud. This com- 
pensation of the spin explains why some of these com- 
pounds do not order magnetically. The main point to be 
explained in the lattice case is the competition between 
the Kondo compensation of the localized spins and the 
magnetic interactions between them. In these materials 
this interaction is mediated by the conduction electrons 
(RKKY-type). Actually, since the Kondo temperature 
is very small it is difhcult to explain why the RKKY 
does not always prevail. Related to this competition is 
the effectiveness of the compensating cloud around each 
/-level. The size of this cloud has been subject of con- 
troversy. Arguments slipw that it should be a large scale 
of the order of vp/Tx^ but other arguments claim to be 
^ a (a is the lattice constant) Jl3 This is a relevant issue 
in the lattice case related to Nozieres exhaustion prob- 
lem which states that there are not enough conduction 
electrons to screen the /-levels. 

To increase the complexity the system may also or- 
der into a superconducting state. Many questions have 
been raised starting from the result that the discontinu- 
ity of the specific heat at is large, of the order of the 
specific heat itself in the normal phase (which originates 
in the heavy fermions). This indicates that pairing oc- 
curs between the heavy /-level electrons, which will then 
form the condensate. Within the Anderson lattice model 
the strong correlations and the hybridization are respon- 
sible for the high effective masses and it has been pro- 
posed that the mechanism for superconductivity lies in 
the strong Coulomb interaction between the /-electrons, 
not in a phonon mediated attraction. 

Using Coleman'sEj slave boson formalism together 
with a large-TV approach, various attempts have been 
made to search for the existence of superconducting in- 
stabilities in |the infinite-t/ Anderson-lattice model. It 
was proposcdtJ that slave bosons fluctuations can pro- 
vide an effective attraction between the electrons to lead- 
ing order in Later, a calculation of the electron- 
electron scattering amplitude to order 1 /iV^ revealed an 
effective attractive interaction in the p and d channels, 
which was interpreted as a manifestation of the RKKY 
interaction, showing that spin fluctuations are an impor- 
tant mechanism.Ej 

Assuming that the normal state is a Fermi liquid, sev- 
eral other studies of superconductivity have been carried 



out-|Cp_the Anderson lattice model and generalizations of 
itoil3 By adding an attractive nearest-neighbor inter- 
action between the /-electrons, so as to explicitly pro- 
vide an attractive channel leading to superconductivity, 
a mean-field study has been carried out as a functiaa. 
of the local repulsion U. Romano, Noce, and Micnas,E£l 
have found a superconducting ground state for finite val- 
ues of U, but no superconductivity was found for large 
values of onsite Coulomb repulsion, in the Anderson lat- 
tice. This is so because the authors consider the Kondo 
regime (this is, ea <^ fj, where /i is the chemical poten- 
tial), where the occupation number of an /—orbital, nj, 
is close to two for small U. Therefore, upon increasing 
the interaction U, this number is reduced to one, blocking 
charge transport in the /-band. 

In this paper we carry out a mean field study of su- 
perconductivity in the U — oo Anderson lattice where 
an attractive interaction between neighboring /-orbitals 
is explicitly introduced in order to simulate an effective 
interaction (which might have various causes) leading to 
superconductivity. Since t/ = oo, we are restricted to 
/-level occupancies in the range < n/ < 1. In the 
mixed valent regime, where n/ is between zero and one, 
charge movement is allowed among the /— orbitals, even 
when U = oo. We study the dependence of the critical 
temperature and /-level ocupancy on the various model 
parameters for different Cooper pairing symmetries. The 
paper is organized as follows: in section || we present the 
model Hamiltonian we use in our study and derive the 
mean field equations. Particular attention is paid on the 
form of the superconducting pairing term. In section III 
we present our calculations of the critical temperature as 
function of the several parameters of the model and we 
summarize our findings in section 



[V 



II. THE MODEL HAMILTONIAN 

We consider an extended version of the Anderson lat- 
tice model, which includes a density-density attraction 
between the electrons occupying neighboring /-orbitals. 
This form of interaction enables us to consider three pos- 
sible symmetries for electron pairing: s, d and p-wave. 
The Hamiltonian is given by 



where 



i,7n 
fc,m 

Hcf=VY,{i,nkm + fl 
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(2) 
(3) 

(4) 
(5) 
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and 



Hj = -J 
^ 2 



E 



(6) 



<i,j> .jn^m' 



where i and j are nearest neighbor sites and rii^™ = 

fimfi,in- The c and / operators are fermionic and obey 
the usual anti-commutation relations. The hybridization 
potential V is assumed to be momentum independent. 
The term Hu represents the strong onsite repulsion be- 
tween the /-orbitals and in the rest of this work we shall 
consider U = oo. The term Hj explicitly describes an 
effective attraction between neighboring /-sites ( J < 0) 
which is responsible for superconductivity. The totaLap=. 
gular momentum projection m takes on N values.E3^E£l 
We shall assume that the local angular momentum of 
the /-sites is half-integer and, therefore, that N is even. 
The term Hj may be re-written in momentum space 

as 



E E 

Q,k,k' 



2 §+fc',m 



(7) 



where the interaction Jr r, 

k,k' 



J ^^exp — k') ■ 6 and 

the summation over S runs over the nearest neighbors. 
Considering the case of a cubic lattice, the interaction 
Jj^ p may be separated into terms with s— , p— and 
d— wave symmetries asS 



boson formalism, we associate a boson operator with ev- 
ery / operator in (|^) in order to prevent double occu- 
pancy at the /-sites. Taking also into account the bo- 
son condensation, we obtain the superconducting part 
of the mean-field Hamiltonian from the substitution: 
/W./ ^ 2;/t/t < zff > +h.c.. FoUowins-these ideas 
we write down the effective Hamiltonian as:Eil 



k,rn 

+ ^^E(4™^-"+4,™4™) 

k.rn 



k.m 



2J 



(9) 



where TVs denotes the number of lattice sites and e f is the 
renormalized energy of the / orbitals due to the on-site 
repulsion. The angular momentum projection m' — —m 
if electron pairing in a singlet state (s- or d-wave) is con- 
sidered and to' = TO in the case of p-wave pairing. The 
gap function ^ = ^fc^m ^^'^ the superconducting or- 
der parameter A^ is given by 



^™ ~ N E ^k ^ f-k,rnfk,rn ^ ' 
k 



(10) 



7 7 1 J^^J"^ 

-^k.k' -■^ % 



where 



where ry^r denotes any of the possible pairing symmetries 
considered in (^). 

The density of the boson condensate z minimizes the 
free energy of the system and £/ is obtained after impos- 
ing local particle (boson-|-fermion) conservation at the 
/-sites: 



- [ cos(fc^) + cos{ky) + cos(fc^) ] , 



is) 
% - 

4^'^' = ^/2 sin(fcO , 

(<i,2_„2) ^ C0S(/C2;) - COs{ky) , 



(d^2_3^2) 



V3 



COs(fca;) -|- COs(fcy) — 2 COs(fcz) 



(8) 



Electron pairing in the superconducting phase will occur 
in the state with total pair momentum Q — 0. 

We implement the condition U =p-po within the slave- 
boson formulation due to Coleman,t3 in which the empty 
/-site is represented by a slave boson bi and the physical 
operator fi in equation is replaced with b^ift- Con- 
densation of the slave-bosons can be described by the 
replacement bi —>< bi >— < fotj ^/z. The mean- 
field treatment of the interaction term Hj involves the 
usual decoupling of destruction and annihilation opera- 
tors but, in keeping with the spirit of Coleman's slave 



z = l- — Y < f- k >, (11) 

k.Tfi 

" '° " '2^,T^ ^lm%^ > + ™ > 

k.m. 



Sea:.-. 



(12) 



Equation (|ll|) states that the mean number of electrons 
at an /-site is 1 — z. 

In order to derive the gap equation and the spectrum of 
elementary excitations we use the Gorkov Green's func- 
tion approach. The anomalous Green's functions that we 
need to consider are: 



^Uk.r- r') = < T./t ^(r)/l- _(r') > , (13) 
^l^f^^-- = < ^Aj-^f-k.-J-'^ > ' (14) 
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and we must also define three other Matsubara Green's 
functions: one that is associated with the conduction 
electrons, another one for the /-electrons and the third 
one is related to the hybridization of the / and c bands: 

Gcmik, t^t') = -< TrC^ Jryl jr') > , (15) 
gf,^{k,r-T') = < Trfj: Jr)fl jT') > , (16) 

GcfMk, T~r')^-< TrCj:jT)fljT') > , (17) 

After Fourier transforming these functions into frequency 
space, we may write down their equations of motion 
(Gorkov's equations) according to the Hamiltonian (||): 

{-iuJn + £/ - fJ-)Gf,rn{k, iuJn) + V^/zQcf,m{k, iuin) 

+ Jz2A„(fc)^)^^(^, -ZC^„) = -1 , (18) 



{-lUJn + Cfe - Ai)^c/,m(fc,«W„) + V^/zgf^rn{k,iuJn) = 0, 

(19) 

{-iljJn ~ e^ + M)-^c/,m(^'*^") ~ ^'^^]^rrXk,iljJn) = 0, 

(20) 



{-iuJn - £/ + f^)^}.mik, i^n) " V y/Ij'lj „^{k , iuJn) 

+ Jz^Al{k)gf^^{-k, zw„) = . (21) 

Diagonalization of the above equations yields the energies 
of the poles of the Green's functions (excitation energies) 
and the corresponding residues (coherence factors). The 
solutions are of the form 



= 1,2 a=± 



iuJn + aEi 



(22) 



The coherence factors, uf and the excitation energies, Ei 
are given in the Appendix. 

The Green's functions have to be determined self- 
consistently using the mean field equations (|ic|)-([l2|). 
These equations can be rewritten in terms of Green's 
functions as 



T 



2 = 1 - ^ E E ^/."'(^' ^^") ^ 



(23) 



' k,m 



£/ - eo = 



VT 



y^y^gc/,m(fc,ia^n) 



k.m 



zJ 
zJT 



(24) 



(25) 



For a given number of particles per site, n, these equa- 
tions must be supplemented with the particle conserva- 
tion condition which yields the chemical potential fi for 
any temperature: 



n = 1 — z 



N, 



(26) 



S 

k'. 



III. RESULTS 

In what follows we consider a cubic lattice in which the 
conduction band dispersion has the simple tight-binding 
form: 



-2t cos(fei 



SO that D = Qt \s half the bandwidth. We have used the 
subroutine hydrd.f from MINPACK in order to solve the 
four coupled equations (p3|)-(p6|). 

The possible pairing symmetries expressed in Eq. (|^) 
have been studied separately. The two rj^: functions cor- 
responding to the cZ-wave symmetry in equation (|^) de- 
scribe different spatial orientations of the angular mo- 
mentum of the Cooper pairs and give degenerate solu- 
tions. This same remark also applies to the three p-wave 
r]^ functions in Eq. (||). 

The critical temperatures Tc are obtained solving the 
mean-field equations using the normal state Green's func- 
tions. On the other hand, the study of A(T), z{T), 
e/(T), and the specific heat requires solving the mean- 
field equations with the full Green's functions. In the nor- 
mal phase, the slave boson condensation temperature, Tz , 
above which z = 0, is given by = (e/ — /z)/ ln(A^ — 1). 
If = 2, z is always finite. For larger values of N , and in 
particular in the limit A^ — > cx), z — > as the temperature 
increases. Correspondingly, — s- 1 and the /-electron 
superconductivity is inhibited. Therefore, for large val- 
ues of A^ it is expected that the mean-field theory will 
not yield superconductivity. One then has to take into 
account the boson fluctuations. We will focus our atten- 
tion in the case N = 2, relevant for instance for Ce and 
Yb materials, but we will return to this point later. 

Figure |l] shows the behavior of the superconducting 
critical temperatures Tc as function of the particle den- 
sity per channel n/N, for each of the three pairing sym- 
metries. It is readily seen that the critical temperatures 
associated with d— and p-wave pairing follow similar 
trends and that the d-wave symmetry exhibits the high- 
est Tc up to densities of about n/N « 0.6. At higher 
densities, a crossover occurs into a regime where the ex- 
tended s-wave pairing becomes the most stable, for the 
parameters considered. 

The value of Tc vanishes at low densities because the /- 
level occupancy also becomes small in that limit (z — > 1) 
and Cooper pairing occurs only between the /-electrons 
in the model under consideration. In the high density 
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limit, Tc vanishes because each /—level is almost fully 
occupied with one electron (z — s- 0), and freezing of the 
charge fluctuations [arising from the term /^/^ in (|^)] 
occurs because of the infinite on-site repulsion. 

Heavy-fermion behavior in the normal phase occurs 
when the chemical potential /i lies close to the peak of the 
density of states (hence the strong effective mass) . This 
peak is the equivalent of the Kondo resonance peak which 
appears in the single-impurity problem. For the lattice 
problem, two strong peaks appear due to hybridization 
between the conduction electron band and the disper- 
sionless band of localized /-states, leading to the large 
electron's effective mass. For densities above n/N « 0.7 
the chemical potential becomes close to the density of 
states peak in the lower band. 

In the superconducting phase the full solution of Eqs. 
(^)-(^) yields a renormalized excitation energy spec- 
trum. In Fig. (H) we show the band structure for 
n/N — 0.7 in the normal and superconducting phases. 
It is clear that for this density /i is in the flat region of 
the band in the normal state. 

It is seen from Fig. |^ that as the density per channel 
n/N approaches 1, the value of Tc is strongly reduced 
until it eventually vanishes. From the same figure one 
can also see that the critical temperature of the s— wave 
state at n/N « 0.7, for instance, is higher than that of 
the d— or p— wave states. For the model parameters con- 
sidered in Fig. |] this means that as the temperature of 
a normal system is lowered, the system would first enter 
a superconducting state with extended s-wave symme- 
try. On lowering further the temperature, the nature of 
the superconducting state becomes a mixture of different 
symmetries. This sequence of phase transitions would be 
different had we chosen different parameters: our calcu- 
lations show that if J/D is less than about 0.4, then the 
critical superconducting temperature of a system with 
n/N « 0.7 would correspond to a d-wave order parame- 
ter (see left panel of Fig. ||). 

The dependence of Tc on the parameters V , eo and J 
shows interesting crossovers. If eo is well below the chem- 
ical potential /i then the /-level is highly populated and 
the system cannot become superconducting unless the 
hybridization parameter V is large enough. On the other 
hand, if eo is not too low a superconducting ground-state 
is obtained even for small values of V . As can be seen 
from the right panel of Fig. |^, Tc first increases with V 
up to a maximum value, but as V is further increased, 
large charge quantum fluctuations at the / orbitals are 
induced and superconductivity is destroyed. Moreover, 
the d— and p— wave superconductivity seem to be more 
stable than the s— wave for large values of V . That this 
result is consistent with Fig. can be easily understood 
as follows: upon increasing the hybridization between the 
/-orbitals and the conduction band, the electron occupa- 
tion in the /-sites is reduced and Fig. |^ already showed 
that depletion of the / band has the effect of reducing Tc 
and increasing the stability of d-wave pairing relative to 
p— and s— wave pairing. 



The temperature dependence of the gap function in the 
superconducting phase is the standard one. In Fig. ^ we 
show a typical case. The crossing of the d— and p— order 
parameters close to Tc is related with the same crossing 
observed in z{T). Since close to Tc, z{T) for the super- 
conducting d— phase becomes slightly higher than for the 
superconducting p— phase, the d— wave phase has an ef- 
fective superconducting coupling that is slightly higher 
than the p— wave coupling, leading to an higher Tc- 

In Fig. 5 we show the dependence of Tc on the /-level 
position. It is seen that the d— wave state has always a 
higher Tc than the p— wave over the range of eo values 
considered. But the s-wave critical temperature exhibits 
a much stronger dependence on cq. In particular, s-wave 
pairing seems to be more strongly depressed for low eo. 

In a normal system at zero temperature the renormal- 
ized /-level energy e/ is located above the chemical po- 
tential and e/ — /i is of the order of the Kondo tempera- 
ture for the equivalent single-imputity problem. Keeping 
the particle density fixed, both ey and ji increase with 
temperature but the difference ef — p, decreases. Our 
calculations show that Tc is smaller than e/ — by a fac- 
tor of about 10 (see Fig ||) over almost the entire range 
of densities considered in Fig. |[ In Fig. || we present 
Tc, e/, and /i for the extended s— wave order parameter 
(the curves for the other symmetries are qualitatively the 
same). The susceptibility dn//de/, in the region of den- 
sities characterized by n/N ~ 0.7 or larger, is very small 
since the /-level density of states is much larger than the 
c-level one, leading to a negative feedback changing the 
chemical potential in such a way as to keep e/ close to 
/i. This is very clear from Fig. |^, where e/ is indeed 
close to for electronic densities where the density of 
states is large. This is consistent with the picture that 
the pairing is developed by the excitations of the system 
resulting from the Kondo compensated lattice. 

Finally we calculate the specific heat for the vari- 
ous symmetries. The non-conventional pairing symme- 
try leads to a power law behavior at low T in the su- 
perconducting phase. In Fig. Q we show the specific 
heat for the various symmetries as a function of tem- 
perature. The specific heat jump at the transition is 
AC/C ~ 1.6,1.3,0.8 for the p, d, s— symmetries respec- 
tively. We have found that the specific heat at low T has 
a T^ dependence for the p, d— symmetries, and has an 
exponential behavior for the s— wave case. 

Considering now the effect of increasing the number of 
channels iV, we find that Tc decreases by one order of 
magnitude or more, as N changes from = 2 to iV = 4. 
For the parameters considered in the figures, the effect 
is most dramatic for s,p— wave symmetries, where su- 
perconductivity is absent for TV > 4. Furthermore, we 
found that the critical temperature of a system with a 
d— wave order parameter is less sensitive to the number 
of channels TV, as compared to the other symmetries. 
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IV. SUMMARY 

Heavy-fermions show a rich and complex behavior at 
low temperatures. In particular, the interplay between 
magnetic correlations, the Kondo effect and supercon- 
ducting correlations is a difficult problem to solve. This 
is further complicated since neither the mechanism nor 
the pairing symmetry are fully established. In this paper 
we have focused on the superconducting order assuming 
that the superconducting correlations are the dominant 
ones. Using a generalized Anderson lattice model with 
nearest-neighbor attraction between the /-electrons and 
with infinite-?/ local Coulomb repulsion, we studied the 
various pairing symmetries using a mean-field approach. 
In this way it is possible to compare the various solu- 
tions in contrast to an approach where, starting from the 
normal phase, the leading instabilities are identified. 

The results show that there are several crossovers be- 
tween the s, d and p-wave pairing symmetries as the pa- 
rameters of the model are varied. In contrast to a previ- 
ous mean-field approach wc find superconducting order, 
even though U — oo. The reason is that we focus on 
a regime where < < 1, while the previous work 
concentrated on a regime where I < nf < 2 (for finite 
U). Since we consider only the case U = oo, Uf has 
to be smaller than one due to the Coulomb repulsion. 
In the previous work as U grows the density n/ ^ 1 
the /-electrons become more localized inhibiting super- 
conductivity. We find the same qualitative behavior as 
we approach the Kondo regime from the mixed valent 
regime. For small values of eo we tend to a regime where 
n/ ^ 1 and superconductivity is suppressed. 

In the mean-field approach ifz — > {uf — *■ 1) the 
gap function A,„ — > 0. This happens for large den- 
sities n/N 1. For Uf ^ superconductivy is su- 
pressed, since the superconducting coupling is among the 
/—electrons. Also, if N is large z ^ at lower temper- 
atures. In particular, p-wave and extended s-wave sym- 
metries are strongly suppressed. For N = 2 z is always 
finite. For larger values of N in general it will be neces- 
sary to consider the boson fluctuations and a treatment 
beyond mean-field will be required. For systems where 
the spin degeneracy is low we expect the results to be 
qualitatively correct. 

We have found that the d-wave and p-wave symme- 
tries yield similar transition temperatures. For large 
nearest-neighbor attraction the extended s-wave pairing 
is preferred. Otherwise, the d-wave symmetry seems to 
be more robust, in particular as N grows. Clearly, we 
are not considering magnetic correlations in our mean- 
field study and therefore the description applies to sys- 
tems where there are no local moments (and therefore 
Tc <Tk) and where Tc > Trkky- 

We found that superconductivity is preferred in a 
mixed valent regime (due to the infinite Coulomb repul- 
sion). There aije-, materials that are mixed valent and 
superconductorsE3EZI. In the framework of weak coupling 



BCS theory one would expect that the local magnetic 
character of the /-states should be pair breaking. How- 
ever, the heavy fermion superconductivity in the Kondo 
limit (integer valent case) reveals that the pairing is of 
another nature that compensates the pair breaking ef- 
fects of the local magnetic character. For materials such 
as CeRu^Si2 there is a considerable mixed valent char- 
acter and accordingly the effective masses are not high. 
Also, the Wilson ratio is close to one indicating a conven- 
tional weak-coupling BCS superconductor. Other mixed 
valent superconductors are not conventional supercon- 
ductors. It would be interesting to identify systems that 
by changing the mixed- valent character could change the 
superconducting temperature, Tc- In the framework of 
our model this would require /-states with large U val- 
ues. The nearest-neighbor attraction could be due to 
several mechanisms like spin fiuctuations or slave boson 
fluctuations (Coulombic nature). 
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APPENDIX: POLES AND COHERENCE 
FUNCTIONS FOR THE GREEN'S FUNCTIONS 

The algebraic solutions of Eqs. ( ^"8[]2l| ) for the Green's 
functions t//,m(fc, ia;„), Qcf.mik^iojn), .Fj ,„(&, iw„), and 
Gc,m{k,iLUn) havc the form 

g{k, ^tOr.) = - E E ■ ' (Al) 

and the coherence factors uf and the excitations energies 
Ei are given below. 

The energies Ei have the form 

E, = y/7/2-/?74^, (A2) 

E2 = ^7/2+7774^ , (A3) 
with 7 and /3 given by 

7 = (£/ - + {ej: + 2V^z + | Jz2A(fc)p , (A4) 
P = [(e,, - M)(e/ - m) - V^zf + I Jz^A(fc)p(,^^ - ^,f . (A5) 

The uf factors for Qf,m{k,iLOn) are given by 



u+ = F{Ei + e^- - m)Ai , (A6) 

u~ = F{Ei - + ^i)Yi , (A7) 

u+ - -G{E2 + ej: - fi)X2 , (A8) 

U2 = -G{E2 - eg + Ai)i"2 , (A9) 
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where the functions Xi and Yi (i = 1, 2) are given by 

+ Ef - zV"^ , (AlO) 



+ Ej - zV^ 

and the functions F and G are given by 
1 1 



F = 



2E,{El-El)' 



G = 



(All) 



(A12) 



2E2{El-El) 

The uf factors for ia;„) are given by 

ut = -Jz''A{k)F{El-{e^-nf), = (A13) 
4 = -J^2A(fc)G(£|-(e£-M)'), 4 = -"2- (A14) 

The factors for Qlf^m{k,iijOn) are given by 

u+ = -V^zFXx , = FYi , (A15) 

u+ = V^GX2 , W2 = -GY2 . (A16) 



The uf factors for Gc,m{k, ioJn) are given by 

ut = FQi , = -FRi , 

U2 = -GQ2 , U2 = GR2 , 

where 



(A17) 
(A18) 



Qi = {e^-li- Ei)\Jz^A{k)\'' + (e/ + Ei)Xi , (A19) 
Ri = {e^-H + Ei)\Jz^A{k)\^ - i-cf + Ei)Yi . (A20) 
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0.03 




FIG. 1. The critical temperature as function of the total 
density of electrons per channel, for the U = oo Anderson 
lattice. The parameters are N = 2, ef = -0.25D, V = 0.2D, 
and J = —0.5D. The hopping integral t = 1 and D — 6t. 
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FIG. 3. Left panel: The critical temperature Tc as func- 
tion of the coupling J. Right panel: The critical temperature 
Tc as function of the hybridization parameter V . The total 
electronic density per channel is n/N — 0.7 and the other 
parameters are the same as in Fig. hi 
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FIG. 2. Band structure of the normal and superconducting 



states, along the direction kx 



at zero temperature. 



= kz in momentum space 
The total electronic den- 
sity per channel is nJN = 0.7 and the other parameters are 
the same as in Fig. |l[ The symmetry of the superconducting 
order parameter is extended s — wa ve. The excitation energies 
El and E2 are given by Eqs. (A2) and (AJ). 
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FIG. 4. Superconducting gap A(r) as functions of the 
temperature for the three symmetries considered. The pa- 
rameters are the same as in Fig. |l| 
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FIG. 6. Critical temperature Tc, renormalized energy level 
€f, and chemical potencial /i as function of n/N. The value 
of e/ — is much larger than Tc. The symmetry of the su- 
perconducting order parameter is extended s— wave, and the 
other symmetries follow the same trends. The parameters are 
the same as in Fig. hi. 
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FIG. 5. The critical temperature Tc as function of the 
/-level bare energy eo- The total electronic density per chan- 
nel is n/N — 0.7 and the other parameters are the same as in 
Fig. I 
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FIG. 7. Specific heat C{T) as function of the temperature 
for the three symmetries considered. The parameters are the 
same as in Fig. ^ 
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